Z2 peak of noise correlations in a quantum spin Hall insulator 



m 

o 

< 
0\ 



I 



a 



I 

o 
o 



> 

in 

^* 

O 
m 



X 



Jonathan M. Edge/'j^] Jian Li,^'^ Pierre Delplace,^ and Markus Biittiker^ 

^ Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 
^Departement de Physique Theorique, Universite de Geneeve, CH-1211 Geneve, Switzerland 

(Dated: April 10, 2013) 

We investigate the current noise correlations at a quantum point contact in a quantum spin Hall 
structure, focusing on the effect of a weak magnetic field in the presence of disorder. For the case of 
two equally biased terminals we discover a robust peak: the noise correlations vanish at B = and 
are negative for B ^ Q. We find that the character of this peak is intimately related to the interplay 
between time reversal symmetry and the helical nature of the edge states and call it the Z2 peak. 

PACS numbers: 72.70.+m, 73.23.-b, 72.10.-d, 85.75.-d 



Measurements of current noise correlations can of- 
fer remarkable new insights beyond conductance mea- 
surements [1]. An example of this is the two-particle 
Aharonov-Bohm effect in which the presence of a flux 
can only be determined by measuring noise correlations 
[51I3]. Quantum spin Hall (QSH) systems, discovered 0] 
after pioneering studies of time reversal invariant band 
insulators [H E], are no exception in this regard. So 
far, various measurements have characterised the QSH 
effect from several aspects. First the quantised conduc- 
tance was measured in a QSH bar [J . Non-local current 
measurements subsequently showed that the currents are 
carried via quantised edge modes [7] . The spin polarisa- 
tion of these edge modes was established very recently 
[5] , thereby vindicating the intuitive picture of the QSH 
effect consisting of two time reversed copies of the quan- 
tum Hall effect. Scanning techniques [51-[TT] have now 
provided additional insights into, e.g., inelastic scatter- 
ing in the QSH systems [T2] . 

To investigate current noise in mesoscopic structures, 
a central element is quantum point contacts (QPCs) 
[131 m]. Theoretical studies of QPCs in QSH systems 
have shown ways to test the properties of the helical edge 
states [HHTT] and determine interaction strengths of the 
edge modes [T5]. Current noise studies have been per- 
formed to distinguish one-and two- particle tunnelling 
processes at the QPC ^^. Correlations between current 
noises have also been investigated to this end |2D], as 
has the effect of interactions on the noise correlations of 
the current which is backscattered from a QPC [21] . One 
question remains open, however, despite its direct experi- 
mental relevance: how do the noise correlations vary with 
a magnetic field that breaks the time reversal symmetry 
(TRS)? In this scenario, the topologically protected edge 
states are singular at zero field; otherwise disorder be- 
comes crucial. This is the question that we address in 
the present letter. 

We investigate the noise correlations in a Hall-bar 
structure with a QPC in a QSH system (see Fig. Jsl. 
Our investigation is based on scattering theory [22] . 
which assumes the edge modes to be approximately non- 
interacting channels. Studies of the edge modes in QSH 
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FIG. 1. Six terminal Hall bar structure with a quantum point 
contact. VP 1 and VP 2 are voltage probes which allow for 
inelastic scattering. 



systems using the Luttinger liquid approach have shown 
this to be a good approximation [T2|. The relevant scat- 
tering matrix in the present setup relating contacts 1 to 
4 is a four-by-four matrix (contacts 5 and 6 will be kept 
grounded). We denote the scattering amplitude for an 
electron coming from contact /3 and going into contact a 
to be Saii ■ In previous work |23j , it has been shown that 
in the presence of TRS: Saa = 0, 5*13 = —5*31, S24 = 
—S42 and Sai3 = Sf3a otherwise. All the off-diagonal 
entries of S are in general non-zero |24) . An immediate 
consequence of the vanishing diagonal entries of S is that 
the equilibrium noise (auto-correlation) at each contact 
is universal - it is proportional to the number of open 
channels connected to the contact [52] , but has no depen- 
dence on the details of the QPC. When TRS is broken, 
the diagonal entries of S become also non-zero, signifying 
the onset of backscattering, and the matrix S is only sub- 
ject to unitarity. Using the scattering theory for coherent 
quantum transport [22| . we can readily write down the 
noise correlations in terms of the scattering matrix. We 
will assume in the following the zero-temperature, zero- 
frequency limit. 

We start with the single-source case, namely we set 
eVi > 0, eV2,3,4 = 0, with V^ the voltage at contact a. 



The cross-correlation noise power is then given by [TJ [55] 
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Pc.p = -j{eVi)\S^i\'\Spi\ 



(1) 



This is the partition noise caused by the sphtting of the 
electronic beam at the QPC and it is non-positive. 

Next we consider the more interesting case of two bi- 
ased contacts. To be specific, we set eVi,2 = eVb > and 
eV3_4 = 0. We will focus on the current cross-correlations 
between the two unbiased contacts (3 and 4). In fact 
the choice of the two contacts to be biased/measured 
is immaterial. In this case, the cross-correlation noise 
power P34 contains not only the partition noise similar 
to Eq. ([1]), but also the exchange noise resulting from 
scattering of two indistinguishable electrons coming from 
two different contacts. It is given by 

P34 = -^(eK))[|^3l|'|^4l|' + 1^321' 1^421' 

+ 5*3]^ 5425*32541 + 53i54253254i . (2) 

The exchange noise, corresponding to the second line 
of the above equation, can carry nontrivial information 
encoded in the phases of the scattering amplitudes and 
manifest it through two-particle interference [25j . This 
distinguishes the exchange noise from other measurable 
quantities to which only scattering probabilities are rele- 
vant, such as conductance and the pure partition noise. 
The total noise power P34 is also negative semi- 
definite, which can be seen by simply rewriting Eq. (pi) 

as P34 = -(e^Fo/Zi) 15^3543 + 5^4544^ Here the uni- 
tarity of the scattering matrix has been used to equate 
1 5|]^ 541 -I- 5^2 542 1 with 1533543 -|- 534 544]. One impor- 
tant implication of the above equation is: in the presence 
of TRS, P34 reaches its maximum (zero) as 533 = 544 = 
[26]; when TRS is broken and backscattering sets in, P34 
generally becomes negative. We call this peak in the cur- 
rent cross-correlations the Z2 peak because it is a peculiar 
phenomenon associated with the form of the scattering 
matrix of time-reversal-invariant topological insulators. 
It is clear that this phenomenon does not depend on the 
choices of biased/measured contacts, since we have made 
no special assumption about the contacts so far. 

Physically, the Z2 peak is a result of an exact can- 
cellation between the partition noise and the exchange 
noise. It is known [U [22j that the partition noise, due to 
the particle nature of electrons, is negative semi-definite, 
whereas the exchange noise, due to the fermionic nature 
of electrons, is positive semi-definite. The two contri- 
butions are not necessarily related in generic cases. Here 
however, TRS and current conservation together demand 
that the two contributions be of equal magnitude. Sim- 
ilar cancellations can occur in two other circumstances. 
In one, both outgoing channels are fully occupied at a 
specific energy. This happens, for example, in a QPC 




FIG. 2. Scattering model for the structure shown in Fig. [5] 
in the presence of a magnetic field and disorder. The setup 
is partitioned so that the overall scattering matrix can be 
obtained from combining the transfer matrices (denoted by 
M). A magnetic field causes a fiux $ in the QPC region and 
non-zero backscattering along each arm of the helical edge 
states. Disorder modelled by random scattering phases as 
indicated in the figure leads to a random distribution of the 
overall scattering matrices. 



based on chiral edge states [571 I2H] where both incom- 
ing channels are fully occupied at the same energy. In 
this case the cancellation is trivial because it merely re- 
flects the absence of current fluctuation in each channel. 
Similarly, a properly-timed mesoscopic two-particle col- 
lider with identical sources can lead to a cancellation of 
the noise correlations [53], as recently demonstrated in 
an electronic on-chip experiment |30| . This case is much 
closer to the present case, in the sense that the currents in 
both outgoing channels are noisy by themselves but their 
correlations vanish identically due to the cancellation. 

Remarkably, the Z2 peak persists even when the in- 
coming channels are subject to strong inelastic scatter- 
ing. To model this scenario we employ two voltage probes 
that are coupled to the two incoming arms, from 1 and 
2, respectively ^31] (see Fig.lsl). For simplicity we assume 
the same coupling strength Tp for the two voltage probes. 
Tp is the probability for electrons in the helical channels 
to enter the additional reservoir connected by a voltage 
probe. The vanishing total net currents in the voltage 
probes require the voltage for both additional reservoirs 
to be 



Vp = 



2 + {2-Tp)T + Tp{l-Tp)T^ 
4 - (TpT)^ 



Vr 



: 



(3) 



^. In the strong couphng limit, 
— T); the cross-correlation P34 



where r=|52i|^ = |5i2 
Tp = 1 and Vp = Fo/(2 
measures coherent contributions from the two voltage 
probe reservoirs instead of the original ones 1 and 2. It is 
clear that the effect of the voltage probes in this limit is 
only to substitute Vq in Eq. (pi) by V^. Such a substitu- 
tion obviously preserves the qualitative structure of the 
Z2 peak. 



Having established that the suppression of the current 
cross-correlation is a robust feature of TRS-preserving 
scattering, we now investigate the effect of a weak mag- 
netic field on the present setup. In real experiments, 
transport measurements on mesoscopic devices normally 
display sample-dependent fluctuations when varying the 
magnetic field, due to disorder [32] ■ In the following we 
will include disorder, and, as a consequence, investigate 
the distribution of the cross-correlation noise power P34. 
We consider two major effects of the magnetic field in a 
generic scenario when disorder is included. The first one 
is the TRS-breaking scattering at the QPC; the second 
one is the backscattering that may occur along each arm 
of the helical edge states before approaching the QPC 
[53]. Here we generally assume that at small magnetic 
fields the lengths of the paths between leads and the QPC 
are smaller than the localisation length of the helical edge 
states [23] such that the current cross-correlation is not 
suppressed simply by localisation. We will model the two 
effects separately, but consistently with scattering theory. 

For the TRS-breaking scattering at the QPC, a min- 
imal model requires an additional loop of helical states 
inserted into the contact area between the two pairs of 
original edge states (see Fig. [2| . This loop is coupled to 
the original edge states in a point-like fashion via TRS 
tunnelling. Electrons encircling the loop accumulate an 
Aharonov-Bohm (AB) phase $ (cf. Ref [23])- To obtain 
a scattering matrix for the combined QPC, we adopt the 
transfer matrix approach (see supplementary materials 
for details). 

The transfer matrix describing the local tunnelling be- 
tween one pair of edge states and the loop states, trans- 
formed from a TRS-preserving scattering matrix, is given 
by 



A/n = 



1 



Vn^ 



UxPcTx tCFxP 



/? = 



1 



\/r^^ 



— s 
r 



(4) 



where t 
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|5'34|, S 
= |^23|. 
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are the conven- 
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tional Pauli matrices. The transfer matrix for the interior 
of the loop reads 



Af <T, = 



7($; 

7*(-$) 



, 7(*) 



e*(*+^'). 



where ^pl and (/?„ are the dynamic phases for the lower 
and upper parts of the loop (see Fig. |2| . We have chosen 
a gauge such that the AB phase only enters the lower part 
of the loop. Transforming the combined transfer matrix 
Mqpc = MqM^Mq, we obtain the scattering matrix for 
the magnetic-flux-dressed QPC: 



S, 



QPC 



P<Jx 

axP 



-Ai($) A2($) 
Ai'(-$) Ai($) 



Pgx 
<JxP 



(6) 



with Ai($) = tax[(Jo + 7(*)7(-*)]/[o-o + t^l{^)-f{-^)] 
and A2($) = (1 - t^)Y^\-<^) -f t2^($)]~i. In order to 
illustrate the effect of the magnetic flux on the scattering 
amplitudes, we extract from Eq. (|6| that 



isin<I> 



533($) = 


' t cos$ + cos((/3-2ilnt) 


534($) = 543(-$) 


= t 


r s2 r2 1 

1 1 1 


" ' f2 + £-*(¥'+*) i2 ^ e-*('^-*) 



(7) 



(8) 



where Lp — ipi + ipu- Clearly, the backscattering at the 
QPC is suppressed when $ = mod tt. On the other 
hand, by using Eq. ^, we find Ps^i^) = -P34(-<I') for 
the present QPC. 

To take into account the backscattering that occurs 
along one arm of the helical edge states between a lead 
and the QPC, we make use of the weak-field-limit result 
obtained in Ref. |23j and write the scattering matrix as 



SssiB) = 



_(-^l-T{B) 



2pii{B) 



r(B) 



r{B) 



^I - r(B)2e-*«(-S) 



(9) 



where t{B) = exp(— ai?^) with a being a sample- 
dependent constant, and £,{B) is a (i3-dependent) ran- 
dom phase. The backscattering can be different for dif- 
ferent arms, which again relies on specific disorder con- 
figurations, but to include all them into the full model is 
straightforward in terms of the transfer matrix approach 
(see Fig. [2]). From this we obtain the full scattering ma- 
trix. 

The full scattering matrix contains parameters that 
are sample-dependent. Varying these parameters allows 
us to obtain distributions of noise correlations as a func- 
tion of magnetic field. For simplicity, we fix t, s, and 
hence r. We also choose a fixed loop area and a fixed a 
that is the same for all arms. The values of these fixed 
parameters are determined as described in the supple- 
mentary materials and they permit a sound comparison 
with numerical simulations that will be presented below. 
At a specific B, we pick randomly the scattering phases, 
namely ipi, tp^ and ^'s, with a uniform probability dis- 
tribution in (0, 27r). This turns out to be sufficient to 
produce a random distribution of the full scattering ma- 
trix (see below). 

The probability distribution of the noise correlation 
P34 produced from the above-described scattering model 
is plotted in the upper panel of Fig. [3] with the overlay- 
ing solid line showing the mean value (P34) as a func- 
tion of B. The Z2 peak of the noise correlation can 
be clearly identified either from the probability distri- 
bution, or more directly in terms of the mean value. 
The peak structure extends from weak magnetic field up 
to the point where the noise correlations are suppressed 
again due to strong backscattering in individual arms. 
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FIG. 3. Probability distributions of the cross-correlation noise 
power P34 as a function of magnetic field B, obtained from the 
scattering model with 100,000 disordered samples in the upper 
panel, and from numerical simulations with 1200 disordered 
samples in the lower panel. Overlaying both panels, the solid 
lines show the mean of P34, clearly showing the Z2 peak at 
B — 0. The dashed lines indicate the average of P34 evaluated 
for a circular unitary ensemble of scattering matrices. 



The maximumly negative value of {P3i{B)) is compared 
with the average of P34 for the circular unitary ensemble 
of four-by-four scattering matrices, given by the dashed 
line. The circular unitary ensemble contains uniformly 
distributed unitary matrices to which the TRS-breaking 
scattering matrices belong [33] . Averaging Eq. (pi) in this 
ensemble yields {P'aa^cue — — (l/15)e'^Vb//i. The maxi- 
mumly negative value of {P^ii{B)) approaches, but does 
not reach, {PaijcuE- This, on the one hand, justifies that 
by only varying the scattering phases a reasonably ran- 
dom distribution of scattering matrices can be obtained. 
On the other hand, it also indicates that this random 
distribution is not quite uniform. 

To examine the validity of our scattering model, we fur- 
ther perform numerical simulations with a microscopic 
Hamiltonian (see supplementary materials for details). 
We construct numerically a device as illustrated in Fig. [5] 
from a lattice model of the HgTe/CdTe quantum wells 3], 



described at low energy by the Bernevig-Hughes- Zhang 
(BHZ) Hamiltonian [6]. Disorder is introduced by adding 
random on-site potentials of a Gaussian profile. The 
resulting potential fluctuation has a magnitude smaller 
than the bulk band gap and a correlation length compa- 
rable to the penetration depth of the edge states. Scat- 
tering matrices connecting transmitting modes between 
leads are computed from Green's functions [34] at various 
magnetic field for each disorder configuration. The prob- 
ability distribution of the noise correlation P34 is then 
obtained by using Eq. (pi), and plotted in the lower panel 
of Fig. I3] Comparing with the upper panel, we observe 
a remarkable agreement between the numerical simula- 
tion and the scattering model at weak field, despite a 
minor quantitative disagreement at stronger field since 



our choice of SssiB) in Eq. (36) is no longer valid. 

To summarise, we have constructed a model for a QPC 
in a QSH system in the presence of disorder, subject to a 
magnetic field. We have computed the ensemble proper- 
ties of the noise correlations for this model and found a 
favourable comparison with results from a numerical cal- 
culation. In particular, both approaches show the pres- 
ence of the Z2 peak, a maximum of the noise correlations 
at zero magnetic field. 
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Supplementary material 



Transfer matrix approach 



We present a derivation of the scattering matrix of figure 2 in the paper (figure H in this supplementary material, 
but with additional labelling) and discuss the number of relevant free parameters. 



General transformation and time reversal symmetry 



Define the scattering matrix S as 



Br - ^ \Ar 



Sll Slr 

Srl Srr 



Al 
ArJ' 



(10) 



and transfer matrix M as 



Ar 
Br 



M 



(Maa Mab 
\Mba Mbb 



Al 
Bl 



(11) 



Here, A and B are vectors of current amplitudes (we will always use A/a for incoming ones and B/b for outgoing 
ones), L and R stand for left and right sides. The two matrices are related by 



M = 
or S = 



'>RL 



' SrrSj^j^SlL SrrSj^j^ 



-MXhMAA 



M 



AB 



Mba - MbbM^]^Maa MbbM^^ 
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(13) 



Mbs, Mo , Mo , Mo M 
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FIG. 4. Partitions of the scattering model. 



In the system we are considering each block S'^,-, M^,- is given by a 2 x 2 matrix. We use the bases defined by 



IhA 




(ai\ 




("A 




(ai\ 


&2 


= S 


0-2 




a^ 


^ M 


a2 


hz 




as 




bs 




bi 


\bij 




\«4/ 




\bij 




\b2) 



(14) 



where ai/bi denotes an outgoing/incoming channel at terminal i. In this basis time reversal symmetry imposes that 
the scattering and transfer matrices satisfy the relations 



S{B) 
M{B) 



"„■ "J *"(-«' To °. 



-a, 
CT. 



M*{-B) 



a, 
-a. 



or 



Srb{B) = -azS'^j^{~B)az, 
Slr{B) = -a,Sl^{-B)a,, 
Maa{B) = a,M*BB{-B)a,, 
Mab{B) = -<7,M*BA{-B)a,., 

where B stands for the magnetic field, L = (1, 2) and R = (3, 4). 



(15) 
(16) 



(17) 
(18) 
(19) 
(20) 
(21) 



Single point contact 



At a single point contact for helical edge states, the scattering matrix is given by 



/6i\ 




/ -t -s r\ fai\ 


62 




-t r* s* 




a2 


&3 




s r* t* 




03 


\bj 




\r -s* t* 0) 




\«4/ 



(22) 



where 1 to 4 label the terminals, and |i|^ + |sp + |rp = 1. We can remove the phases of the scattering amplitudes by 
redefining the current amplitudes (which respects time reversal symmetry) as follows 



63 



/^tsrbl\ 
V ^rb4 J 



0-2 _ ^srO.2 
^3 irtO-i 



(23) 



where (,x '■= x/\x\. Then the scattering matrix becomes 



S = 






Since we have removed the phases, r, s,t are now real, non-negative numbers. To test (time reversal) symmetry 
relations, it will be useful to notice that <Jx<yzP<^x<^z = li- 
lt is easy to find the transfer matrix, which will be denoted by Mq (see Figure [i]), for the above scattering matrix 
to be 



&3 



tf3<Jx /3 
Vl - t^ \crxl3ax taxP 



02 
61 

\b2j 



(25) 



Point contact including a loop 



The transfer matrix for the loop, which will be denoted by Af$ (see Figure p|, is given by 

/«3'\ 



0-2" 



7($; 

7*(-$) 



b,' 



7(*) = » ,.^„ 







(26) 



where (pi and (/Ju are the dynamic phases for the lower and upper parts of the loop, $ is the AB phase. We have 
chosen the gauge such that the AB phase only enters the lower part of the loop. Note that the phases appearing 



in Eq. (23) and immediately connected to the loop can indeed be absorbed into </?; and (/?„, therefore the removal of 



these phases in Sec. [] is justified. 

The transfer matrix for the full point contact including the loop is given by 



= Mqpc 



vvy 



(aA 



"2 



= MqA'UMq 

1 //3a, \ /Ai($) A*(-$)\ f/3ax 
l-t^\ cTxiSj ^^A2($) A*i(-$)y' ^^ a,/3 

Ai($) = to,[7^($) + 7*(-*)], A2($) = 7^($) + t'7*(-^) 



M, 



QPC 



The corresponding scattering matrix is given by 

Sqpc 



Pox 
a,/3 



Ai($) = Af(-<i>) = 



A2($) = (l-f 




-Ai($) A2($) 
Ai^(-$) Ai($) 



^/3cr:r 

(7,/3 





t(l+e''^+) 



t(l+e" 







,l+t2e'*+ 
= $ - (/3; - (^„, 0+ = $ + (/3; 







'/'u 



(27) 

(28) 
(29) 

(30) 
(31) 

(32) 
(33) 



For the lower-right block of Sgpc, we find 



^33($) = ^44($) 



rs 



i sin$ 



t cos $ + cos{ipi + tpu ~ 21 In t) 



^34(*) = ^43(-$)-t 



1 



^2 _|_ g-i(ip,+i^u+*) ^2 ^ g-i((^,+ipu-*) 



(34) 

(35) 



These are the same as equations (7) and (8) in the paper. In a similar way, the other entries of the scattering matrix 
can be obtained. 



Including back-scattering 

Unitarity |33j constrains the scattering matrix for back-scattering, using the arm at contact 1 as example, to 






(36) 



where p, t, ^i and ^2 are in general all functions of magnetic field _B; r > 0, pp* + r^ = 1, |^i| = |^2| = 1- The 
principle of micro-reversibility imposes further constraints: p{—B) = —p{B), t{B) = t{—B) and ^2{B) = ^i{—B). In 
the following we write p = \p\^p with ^p{—B) = ~^p{B). 



The transfer matrix corresponding to Eq. ( 36 1 is 



1 flpl^pCi Ci 



6 
p\Cp 



\p\CA 



1 



6 
Q 



1 fCi o\ f\p\^p 1 



6 



1 



\p\C 



ip^^pCiC2,U-B) = -UB)- 



(37) 

(38) 

(39) 
(40) 



The above expressions allow us to remove phases ^1 and ^2 because changing phases of the final current amplitudes (i.e. 
fli's and 6i's) will not lead to any physical effect. Note also that the phases appearing in Eq. (23) and immediately 
connected to the back-scattering part can be absorbed into £^p, therefore the removal of these phases in Sec. I] is 
justified. 

The back-scattering transfer matrices for the whole setup are then in general 
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BS 
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\b2j 






/PI Q 1. 
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T3 

VO 
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\h2J 






Tl 
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•14 
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•14 




■'4 



'A h"\ 






(41) 



(42) 



where r^'s are all real positive and p^'s are complex; pi{—B) — —pi{B), Ti{B) = Ti{—B). 




FIG. 5. Numerical simulation setup. The total size of the setup is l/imx0.6/im. The separation of the quantum point contact 
is 80nm. 



Assuming Ti{B) ~ t{B) ~ exp{—aB^) for all i, we have 



M 



(1.2) 
BS 



1 



P(l,2) CTQ 



M 



r \ CTo P(i 2) 

(3,4) _ 1 fP{3A) CTQ 
BS 



T \ CTq P(3,4), 

P(i.2)(i?) = -P(i.2)(-B) = yr^^ 1^%' l^) , 
P(3,4)(i?) = -P(3.4)(-^) = v/^^^ (^;;' ^^J • 



Combining everything, we have 



Mtotai = M^ig^^MQPcM'i'^\ 



(43) 
(44) 
(45) 
(46) 

(47) 



From this the total scattering matrix can be obtained using equation (13). In total, we have 6 independent random 
phase factors: e*"^' , e"^" and fp^ («=!, 2, 3, 4); the first two are i?-independent, the last four satisfy <^p. {—B) — —£,p- [B). 



Numerical simulations 



In our numerical simulations, we deal with a setup as shown in Fig. [5] The colors in the interior of Fig.lSlencode the 
electric potential fluctuation caused by disorder. Here, only one speciflc example of disorder configurations is shown. 



Hamiltonian 



In the clean limit, we take the Bernevig- Hughes- Zhang (BHZ) Hamiltonian with a spin-orbit coupling term owing 
to bulk inversion asymmetry [H [SJ [SSI • This Hamiltonian reads 



H{k) = C - Dk^ -f {M - Bk^)T^ (?) (Tq + Ak^r.^ ^a^^ AkyTy ® a^ + Ary ay , 



(48) 



where r's and cr's are Pauli matrices corresponding to orbit and spin degrees of freedom respectively. The parameters 
are taken from Ref. [35]. The orbital eff'ect of the magnetic field is taken into account by substituting k with 
k + {e/h)A. where A. is the vector potential. In numerical simulations, we discretized this Hamiltonian with a lattice 
spacing oq = 5nm. 
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FIG. 6. Local density of states as a consequence of biased contact 1. 



Disorder potential 



The disorder potential, for each impurity center r,;, is modeled by a Gaussian function: 

U^{r) = u,exp[-(r-r,)VA^], 



(49) 



where Ui is randomly taken in a uniform distribution in (—Uo/2, C/o/2), A stands for the range of the potential and is 
fixed for all impurity centers. We also introduce a parameter p for the density of impurity centers. The correlation 
function for the resultant potential fluctuation is 

(f/(r)C/(r + 6r)) - ^U^{p\') exp{-6ry2\^) (50) 

where U{r) = J^i Ui{r). In our simulations, we take Uq = lOmcV, A = 25nm, and p\^ = 0.25. 

Scattering matrix 

In each normal metal lead connected to the quantum spin Hall insulator sample, there exist a number of transmission 
modes to ensure a good contact between the lead and the sample. The scattering matrix relates the current amplitudes 
for these transmission modes. It can be obtained from the Green's functions t34i,i36i: 



S = t- i^/AW'^G^W^ 



(51) 



where G^ = {Ef + iri-H, 



sample " 



-WGf^^^gW^)'^ is the retarded Green's function for the sample, Gf^^^^ is the retarded 
Green's function for the semi-infinite leads in their eigen-mode basis, A = —2Iin{Gl^^^^^) is the spectral function which 
is diagonal and positive semi-definite, and W is the coupling between the eigen-modes of the leads and the sample. 
In our simulations we take Ef — 0.5meV, where the edge state penetration depth is about 50nm. 



Comparison with scattering model 



In order to make a sound comparison between the numerical results and the scattering model based on edge states, 
we determine the fix parameters used in the scattering model in the following way. The area of the loop is estimated 
from a typical picture of the local density of states around the quantum point contact of our setup (see Fig. ro| . It 
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takes the value 0.012/ini^. t, s and r are chosen such that, in the absence of magnetic field, the overall scattering 
probabilities for the QPC obtained by averaging the scattering phases ipi and (pu equal those obtained from numerical 
simulations. Finally, the factor a is evaluated from a typical dependence of the total transmission probability (from 
one terminal) on the magnetic field. It takes the value 6T^^. 

I 



